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Although structures in the Universe form on a wide variety of scales, from small dwarf galaxies to large super clusters, the 
generation of angular momentum across these scales is poorly understood. Here we investigate the possibility that filaments 
of galaxies—cylindrical tendrils of matter hundreds of millions of light years across—are themselves spinning. By stacking 
thousands of filaments together and examining the velocity of galaxies perpendicular to the filament's axis (via their redshift 
and blueshift), we find that these objects too display vortical motion consistent with rotation, making them the largest objects 
known to have angular momentum. The strength of the rotation signal is directly dependent on the viewing angle and the 
dynamical state of the filament. Filament rotation is more clearly detected when viewed edge-on. In addition, the more massive 
the haloes that sit at either end of the filaments, the more rotation is detected. These results signify that angular momentum 


can be generated on unexpectedly large scales. 


text is one of the key unsolved problems of cosmology. In 

the standard model of structure formation, small overdensi- 
ties present in the early Universe grow via gravitational instability 
as matter flows from under- to overdense regions. Such a potential 
flow is irrotational or curl free: there is no primordial rotation in 
the early Universe and angular momentum must be generated as 
structures form. 

Tidal torque theory'” provides one explanation—the misalign- 
ment of the inertia tensor of a gravitationally collapsing region of 
space with the tidal (shear) field can give rise to torques that spin 
up the collapsing material’**. Such an explanation is valid only in 
the linear regime, namely in the limit where density perturbations 
are small with respect to the mean and where flows are laminar. 
As a collapsing region reaches turnaround, tidal torques cease to be 
effective and the final angular momentum of a collapsed region is 
far from what tidal torque theory would predict”. Although one 
recent study’ has demonstrated that galaxy spin direction (that is, 
clockwise versus anticlockwise) can be predicted from initial condi- 
tions, revealing a critical clue to the nonlinear acquisition of angular 
momentum, our understanding of spin magnitude, direction and 
history remains in its infancy. Regions that are still in the linear or 
quasilinear phase of collapse could provide a better stage for the 
application of tidal torque theory. 

Cosmic filaments’, being quasilinear extended topographical 
features of the galaxy distribution, provide such an environment. 
Yet, owing to the challenges in characterizing and identifying such 
objects, potential rotation on the scales of cosmic filaments has been 
discussed“ but never measured until now. 

It is known that the cosmic web in general and filaments, in 
particular, are intimately connected with galaxy formation and 
evolution'®'*. They also have a strong effect on galaxy spin'””, 
often regulating the direction of how galaxies**’ and their dark 
matter halos rotate**-*!. However, it is not known whether the cur- 
rent understanding of structure formation predicts that filaments 
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themselves, being uncollapsed quasilinear objects should spin. A 
recent study (published on arXiv while this draft was being final- 
ized)” examined the velocity field around galactic filaments defined 
by halo pairs in a large N-body simulation and found a statisti- 
cally significant rotation signal. This is an intriguing finding and, 
although filaments and their rotation speed are defined differently, 
the current work in which the observed galaxy distribution is exam- 
ined in a bid to find possible filament rotation was partly motivated 
by the theoretical suggestion that filaments may spin”. 


Results 

After segmenting the galaxy distribution into filaments using a 
marked point process known as the Bisous model”, each filament 
can be approximated by a rectangle on the sky and thus the galaxies 
within it may be divided into two regions (A and B) on either side 
of the filament spine. The mean redshift difference Az,, of galaxies 
between two regions are considered as a proxy for the line-of-sight 
velocity difference and hence for the filament spin signal. Since 
measuring mean values among subsamples will by construction 
result in differences between such averages, any measured value of 
AZ, needs to be assigned a significance based on a randomization 
procedure (explained in Methods). Figure 1 shows the statistical 
significance of the measured Aza, as a function of z,,,,,/AZ,3 (where 
Z,ms, is the root mean square (r.m.s.) of the galaxy redshift)—a proxy 
for the dynamical ‘temperature of the filament (Methods). The 
number of galaxies in each region is denoted by colour (the effect 
of galaxy number on this signal is shown in Supplementary Fig. 1). 
Two salient points can be gleaned here. First, the more galaxies in 
a given filament, the more inconsistent the redshift difference Az,, 
is with random. Second, (as expected) the colder the filament, the 
more inconsistent the redshift difference is with random. This sec- 
ond point is a generalization—cold filaments with z,,,,./AZ,,3<1 
and few galaxies can have redshift differences only weakly incon- 
sistent with random expectations. However one may note that as 
a trend, the colder the filament, the more significant the redshift 
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Fig. 1 | The statistical significance of filament rotation. The statistical 
significance of Az,, being consistent with random is shown as a function 

of the filament dynamical ‘temperature’, z,,,./AZap, in which AZ, is the 
redshift difference of galaxies between the approaching and receding 
regions of each filament. The higher this quantity, the more unlikely it is that 
Az, is a random occurrence. Each filament is coloured-coded by the number 
of galaxies it contains, with yellow being rich filaments and purple being 
poor filaments as designated by the colour bar. At a given value of Z,ms/ 
AZ,g the AZs seen in richer filaments is more statistically inconsistent with 
random. ‘Cold’ or ‘hot’ filaments are separated by z,,,,/AZ,,= 1. Note that 


No is the galaxy number in either region A or region B. 


difference. In other words, if Az,, is considered a proxy for filament 
spin, we observe a spectrum of filaments from dynamically hot that 
are consistent with random to dynamically cold filaments that are 
completely inconsistent with random at the many-sigma level. Note 
that even for dynamically hot filaments, there are a few that are 
highly inconsistent with random. 

The reader will note that the wide distribution of significance 
seen at a given value of z,,,,,,/AZ,p, is also a reflection of the distribu- 
tion of inclination angle made by the filament axis with the line of 
sight. The median value of z,,,,/AZ,, (as well the standard devia- 
tion) is shown in Fig. 2 as a function of the angle subtended between 
the filament spine and the line of sight. Two important points can 
be observed by this plot: (1) the filament’s dynamical temperature 
(Zms/AZap) increases as filaments become parallel to the line of 
sight. Because the actual value of z,,,,/AZ,, for a given filament 
is independent of the inclination angle ¢, this fact demonstrates 
that the viewing angle is critical in measuring this quantity and that 
Zrms/AZaz can indeed be used to measure rotation for filaments 
that are close to perpendicular to the line of sight. Furthermore 
such a relationship is natural if a significant component of a galaxy’s 
velocity is perpendicular to the filament spine and consistent with 
rotation or shear. (2) The standard deviation of Z, ms/AZ,g increases 
as filaments become parallel to the line of sight. This is expected if 
in addition to rotational motion about the filament spine, galaxies 
are also travelling along the filament axis—that is, helical motion. 

The cumulative distribution of Az, for both the entire observed 
filament sample, as well as the randomized trials, is shown in 
Fig. 3. The cumulative distribution is shown for all filaments 
(Fig. 3a), filaments whose axis is inclined by cos <0.5 to the line 
of sight (Fig. 3b) and dynamically cold filaments (Z,,,,,/AZap< 1) 
inclined by cos# <0.2 (Fig. 3c). The reader will note that even when 
examining all filaments where the inclination angle is completely 
ignored (and hence includes filaments viewed along their axis that 
will probably weaken the signal), the full distribution of Az,,is 
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inconsistent with randomization tests (see Supplementary Fig. 1 
for an examination of the effect of filament inclination angle). This 
inconsistency increases when considering filaments with cos < 0.5 
and cos <0.2 that are cold. To quantify the statistical significance 
of the cumulative distribution one may simply measure, in units of 
the 10,000 random trial’s standard deviation, how far the measured 
signal is from the mean randomized signal. This is plotted in the 
top panels of Fig. 3 and shows that randomized trails are statisti- 
cally inconsistent with the measured signal at very high confidence. 
In other words, shuffling the redshifts of galaxies in a filament is 
unlikely to produce redshift differences as great as that observed. 

The redshift difference (a proxy for the rotation signal) stacked 
across various (sub-)samples, is presented in Fig. 4. The following 
conventions have been adopted. Region A (defined as the region 
with greater mean redshift) is plotted in the top part of each plot and 
region B is plotted in the bottom part. The position of each galaxy 
along and perpendicular to the filament axis is shown on the x and y 
axis, respectively, and along the x axis is normalized to the filament’s 
length. Each galaxy is coloured by its redshift difference Az, with 
respect to the mean redshift of all galaxies in the filament according 
to the colour bar on the right. In the ideal situation where all galax- 
ies exhibit circular or helical motion about the filament axis, such 
a plot would only have red points in the top part and blue points in 
the bottom part. The statistical significance of each (sub-)sample, in 
units of ø, is indicated on top of each panel. 

In Fig. 4a,d,g, we show the stacked rotation signal for all fila- 
ments, filaments whose axis subtends an angle cos <0.5 with the 
line of sight and filaments whose axis subtends an angle cos < 0.2 
with the line of sight, and that have z,,,,/Az,;<1. Comparing 
Fig. 4a with Fig. 4d, the reader will note what has been mentioned 
before, namely that merely changing the inclination angle increases 
the signal. Figure 4g shows a very strong rotation signal at 3.30— 
when considering dynamically cold filaments that are mostly 
perpendicular to the line of sight, the rotation signal becomes 
very convincing. 

Since filaments are long tendrils of galaxies often connecting 
nodes of the cosmic web, the mass of the two halos closest to the 
filament’s two end points are examined for a possible correlation 


Angle ¢ to line of sight (°) 
90 75 60 45 30 0 


3.0 4 


2.5 4 


2.0 7 


Z ms/^ ZiB 


cos(#) 


Fig. 2 | The median dynamical ‘temperature’, z,,,,,/AZag, of filaments as 

a function of the inclination angle p between the filament spine and the 
line of sight. z,,,,. is the r.m.s. of galaxy redshift and Aza, is the redshift 
difference of galaxies on either side of the filament spine. That this quantity 
increases as the filaments become parallel to the line of sight is consistent 
with galaxies moving predominantly around the filament spine. Error bars 
represent the standard deviation about the median. 
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Fig. 3 | The cumulative distribution of Az,,, the redshift difference of galaxies in the approaching and receding regions of each filament. a, All filaments 


irrespective of viewing angle. b, Filaments whose axis subtends an angle cos < 0.5 with the line of sight. c, Filaments whose axis subtends an angle 


cosġ < 0.2 with the line of sight and that are dynamically cold, namely z 


r.m.s. 


/ AZ, <1. The red solid line shows the distribution of observed filaments and the 


10,000 grey lines indicate the distribution after shuffling the redshifts of galaxies randomly. The median value of these 10,000 random samples is shown as 
the black solid line. The top panels measure, as a function of Az, the distance, in units of the standard deviation of the randomized distribution, between the 
measured curve and the mean of the random distributions. The top x axis displays the rotation velocity of the filament in kms“, calculated as y =c X AZyp. 


with the signal strength we measure. We thus perform the following 
procedure. We search the group catalogue“ for the two groups that 
are closest to the filament’s two end points. The reader is referred to 
that paper for details. In sum, the galaxy groups are detected with 
a modified friend-of-friend algorithm, where the friend-of-friend 
groups are divided into subgroups so that each subgroup can be 
considered as a virialized structure. The mass of each group is esti- 
mated using the virial theorem, the size of the group in the sky and 
the group velocity dispersion along the line of sight. By summing 
these two group masses, we can ascribe a single mass abutting each 
filament. We call this mass the ‘filament end point mass. 

The three samples of filaments can then be subdivided according 
to this characteristic. The middle and right columns of Fig. 4 show 
the stacked signal for the 10% of the filaments with the smallest fila- 
ment end point mass (that is, Fig. 4b,e,h) and largest 10% filament 
end point mass (Fig. 4c,f,i) associated halo end points, respectively. 
The limits of the 10th and 90th percentile group masses of Fig. 4h,i 
are ~10'°4 and ~10'*%, in units of Mo, respectively. 

A clear dependence on filament end point mass is seen for all 
three samples in Fig. 4. The smaller the filament end point mass, the 
weaker the rotation signal. The filaments with the largest end point 
mass, most perpendicular to the line of sight, and dynamically cold, 
show an outstanding 4.20 inconsistency with the null hypothesis of 
random origin. The fact that the mass of the clusters at the end of 
each filament directly affects the measured signal is evidence that 
the signal we are measuring is a physical effect related to gravita- 
tional dynamics. 

A ‘rotation curve’ can be constructed for the stacked filament 
sample by converting the redshift difference into a velocity dif- 
ference computed as cx Az, where Az is the redshift difference 
between all the galaxies at given distance and the mean redshift of 
all galaxies in the filament, and c is the speed of light. This is shown 
in Fig. 5. We adopt the typical convention that negative (positive) 
velocities corresponds to approaching (receding) motion. Figure 5 


shows that rotation speed increases with distance, reaching a peak 
of around ~100kms"! at a distance of ~1 Mpc. It then decreases 
with distance, tending to zero at a distance of greater than ~2 Mpc. 
Incidentally, this behaviour motivates the choice of filament thick- 
ness in this study, and should be kept in mind in future studies of 
the dynamics of cosmic filaments. 


Summary and discussion 

The results presented in this paper are consistent with the detection 
of a signal one would expect if filaments rotated. What is measured 
and presented here is the redshift difference between two regions on 
either side of a hypothesized spin axis that is coincident with the fil- 
ament spine. The full distribution of this quantity (Fig. 3) is incon- 
sistent with random regardless of the viewing angle formed with the 
line of sight but does strengthen as more perpendicular filaments 
are examined. It also strengthens when considering dynamically 
cold filaments—filaments whose galaxies have small r.m.s. values of 
their redshift. We emphasize that this is not a trivial finding. There 
is no reason—besides an inherent rotation—that filaments that are 
dynamically cold should exhibit such a signal. 

One may draw an analogy with galaxies here: a disk galaxy 
viewed edge-on shows a strong Az,, (consider Az,, as the Doppler 
shift from stars on the receding and approaching edge of the gal- 
axy). The value (and significance as measured by the likelihood that 
it is random) of such a Az, would decrease as the viewing angle of 
the galaxy is increased until the galaxy is seen face on. At this point 
the galaxy would appear ‘hot’—the motion of stars along the line 
of sight being random. In addition, there is a degeneracy between 
morphology and shape in such situations: a very thin galaxy viewed 
with a large inclination angle appears ellipsoidal. Since not all galax- 
ies are rotationally supported, some galaxies that appear ellipsoidal 
will be inclined rotating thin disks while others will be dispersion 
supported ellipsoids with little inherent rotation. It is this combina- 
tion of inclination angle and inherent rotation which explains the 
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Fig. 4 | The stacked rotation signal of filaments. a-i, Galaxies’ position r,,, along the filament are scaled by the length r, of the filament. d,,, is the distance 


of galaxies to the filament axis. Region A (defined as the region with greater mean redshift) is shown in the top part of each plot and region B is shown in 


the bottom part. The rows show the stacked rotation signal for all filaments ( 
(d-f), and filaments with cos < 0.2 and that have z,,,;/AZag <1 (g-i). The fil 


rms, 


a-c), filaments whose spine subtends an angle cos < 0.5 with the line of sight 
ament sample is divided up according to the group mass at the filaments’ two 


end points: b,e and h (¢,f, and i) show the stacked rotation signal for the 10% least (largest) filament end point group mass and a,d and g show the signal 


irrespective of this quantity. The redshift difference is indicated in the colour 


spectrum of Z,,,/AZ,p and significance seen in both the galaxy dis- 
tribution as well as the spin of filaments. 

We note that it is the unique combination of small z,,,, and 
large Az,, that returns a strong statistically significant signal. In 
other words, there are filaments with large z,,,, and small Az, 
that clearly do not show any signal. This work does not predict that 
every single filament in the Universe is rotating, rather that there 
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bar. 


are subsamples—intimately connected to the viewing angle end 
point mass—that show a clear signal consistent with rotation. This 
is the main finding of this work. Again, we note that it is not trivial 
that filaments with low values of z,,,, exhibit large values of Az,,. 
This is only expected if filaments rotate. In Supplementary Fig. 2, 
we considered other possible biases that could systematically mimic 


the signal found here, such as the effect of the large-scale motion 
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Fig. 5 | The rotation curve of filaments. The filament rotation speed as 

a function of the distance between galaxies and the filament spine. The 
rotation speed is calculated by cx Az, where Az is the redshift difference 
of galaxies at given distance with respect to the redshift of the filament. 
The distance of galaxies from the filament spine in the receding region is 
displayed in red and ascribed positive values, while the distance of galaxies 
in the approaching region is marked in blue and ascribed negative values. 
Error bars represent the standard deviation about the mean. 


of galaxies, finding these are unable to account for the detection 
found here. 

When divided into subsamples based on the mass of the halos 
that sit at the ends of the filaments, those filaments with the most 
massive clusters have more significantly spinning signals than either 
the full sample or the filaments abutting the least massive halos. This 
gives a hint towards what may be driving the filament spin signal. 
It is also consistent with the idea that the spin is due to tidal fields 
influenced by the presence of massive aspherical density perturba- 
tions. The existence of a mass dependency is also strong evidence 
that what we have observed is not a systematic effect due to (for 
example) cosmological expansion (Supplementary Fig. 2). Any such 
systematic effects on the measurement of the redshift difference and 
its association with a velocity difference would not manifest as a 
mass dependency. The mass dependency opens the door to a num- 
ber of important questions as to the origin of this rotation or shear, 
questions that will be examined in future (theoretical) work such 
as: how is the measured effect engendered by the formation of clus- 
ters abutting filaments? Is there an angular momentum correlation 
between spinning filaments and nearby clusters? And perhaps most 
importantly: at which stage in structure formation do filaments spin 
up, and how does this affect the galaxies within them? 

According to the Zel’Dovich approximation”, the motion of gal- 
axies is accompanied by larger-scale matter flows**’. The mass and 
velocity flows are usually transported in a ‘hierarchical’ way: at first, 
matter collapses along the principal axis of compression forming 
great cosmic walls. Matter then flows in the wall plane along the 
direction of the intermediate axis of compression to form filaments. 
The final stage is the full three-dimensional collapse of an aspheri- 
cal anisotropic density perturbation wherein matter collapses and 
flows along the filament axes to form clusters. Under such a model 
of mass flow, one may expect that filament spin is formed at the 
second stage from the compression along the intermediate axis of 
compression. In other words, we expect motion along the filament 
to still be linear while motion perpendicular to the filament spine 
to be nonlinear. 

The so-called quasilinear regime describes the scales we are 
examining, where the linear approximation is no longer fully able 
to describe the dynamics, but the system is not yet fully nonlinear. 
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‘Shell crossing’ is generally considered to be the ‘moment’ where the 
linear (laminar flow) approximation ceases to accurately describe 
the dynamics, and nonlinear effects start to dominate. In objects 
such as the filaments examined here, shell crossing has occurred 
along the two axes perpendicular to the filament spine but not 
yet along the filament axis itself. Accordingly, a simple analytical 
calculation suggests that when a (significant) velocity difference is 
observed along an axis that has shell-crossed (that is, for filaments 
that are aligned perpendicular to the line of sight), what is measured 
is a vortical flow, and filaments can be said to ‘spin. The term ‘spir is 
a simplification but not an incorrect one. Because a fluid is modelled 
as a vector field, it is customary to use the vorticity to describe any 
rotational motion. When filaments are viewed along the axis that 
has not yet shell-crossed (that is, for filaments that are aligned along 
the line of sight) any motion is ascribed to a shear that, together 
with the vorticity, leads to helical motion. In other words, the signal 
we measure is both consistent with a shear caused by galaxy motion 
in a sheet, as well as the vortical spin of a filament embedded in that 
sheet. These two views are not contradictory with each other, since 
in the picture described above the helicoidal spin of a filament is 
engendered by the shear of the sheet in which it is embedded com- 
bined with vorticity in the directions that have shell-crossed. 

We know that from studies of the spatial distributions of smaller 
galaxies that inhabit the regions between large galaxies that they are 
accreted in a twofold process: first onto the filament connecting the 
two larger halos and then second an accretion along the line con- 
necting the halos. In other words, galaxies are first accreted from 
a direction perpendicular to the filament spine“ and then along 
the spine**°"°. Assuming filaments can be approximated by a cylin- 
drical potential, such a flow would naturally give rise to an rota- 
tion about the filament axis, followed by helical accretion. This is 
consistent with the theoretical picture formed by ref. °, in which 
the analysis of numerical simulations indicated that the angular 
momentum of filaments is generated by those particles that collapse 
along the filament’s potential. Taken together, the current study and 
ref. ® demonstrate that angular momentum can be generated on 
unprecedented scales, opening the door to a new understanding of 
cosmic spin. 


Methods 


The method presented here is fairly straightforward. We begin by using a publicly 
available catalogue of filaments constructed from the Sloan Digital Sky Survey 
Data Release 12 galaxy survey”'*’. Such a catalogue is built using the Bisous 
algorithm”, which identifies curvilinear structures in the galaxy distribution. 

The Bisous algorithm is not astronomy specific and was developed to identify 
roads in satellite imaging**-*’. Each filament is approximated by a cylinder with 

a given axis that makes a viewing angle # with the line of sight. Our hypothesis 
rests on the assumption that if filaments rotate, then there should be a statistically 
significant component of a galaxy’s velocity that is perpendicular to the filament 
spine. Furthermore, this component should have opposite signs on either side of 
the filament spine: one receding and one approaching. To detect such a signal in 
principle, we would only want to examine filaments whose spine is perpendicular 
to the line of sight such that a galaxy’s redshift (a proxy for the radial component 
of its velocity) would coincide with its angular velocity. However, restricting the 
filament sample this way dramatically reduces the sample size. Therefore, we divide 
the filament catalogue according to viewing angle into three bins: all filaments 
irrespective of viewing angle, those inclined by more than 60° (cos# < 0.5) with 
respect to the line of sight, and those dynamically cold (see below) and with a 
viewing angle more than ~80° (cos# < 0.2). This last sample is chosen for the 
reason that filaments inclined by more than ~80° are perpendicular enough that 
the measured value of Z, m.s./AZas is probably a good approximation to the inherent 
one. A consequence of this is that there are more cold filaments with Z,ms/AZ4g < 1 
than hot filaments with z,,,,/AZ,,> 1 inclined by cos# < 0.2. In Supplementary 
Section 1, the effect of viewing angle is examined in more detail. We find that the 
filaments whose spines are closer to perpendicular to the line of sight have more 
significant rotation signals. 

Each filament axis delineates two regions on either side of the filament spine. 
The Bisous filament’s have an inherent scale of around 1 Mpc. However, all galaxies 
that are within a 2 Mpc distance from the filament axis are considered to be in one 
of these two regions. Galaxies farther afield are not considered. We note that most 
of the galaxies that are within 2 Mpc are within 1 Mpc and that including galaxies 


843 


ARTICLES NATURE ASTRONOMY 


farther away from the filament spine will ‘dilute’ the signal as the regions farther 
from the filament axis will, by construction, include random interlopers. 

The mean redshift of galaxies in both regions is calculated. The region with 
the greater mean redshift is arbitrarily called region A and the region with smaller 
mean redshift is called region B. The redshift difference AzZ,,= (Z,) — (Zp) is 
computed. The relative speed is then simply y =cAz,,, where c is the speed of 
light. Next we examine whether the filament is dynamically hot or cold by using 
the redshift z of all galaxies in each filament to compute the r.m.s. namely Z, m.s. 
This is then compared with Az,. Dynamically hot filaments have z,,,, > Az, while 
cold filaments have r.m.s. values less than the redshift difference. Such a definition 
precludes the ability to search for a signal consistent with rotation in hot filaments 
since it is by definition flawed to ascribe a physical significance to a Aza, if it is less 
than z,,,,,- Thus, we do not expect to be able to find a discernible rotation signal in 
dynamically hot filaments. 

A lower limit on the number of galaxies in each region of A and B of each 
filament must be assumed. Without such a limit, we would introduce pathological 
cases into the analysis, namely filaments with no or just one or two galaxies in 
region A or B. Therefore, we only consider filaments with at least three galaxies in 
each region—it is difficult to justify the computation of an r.m.s. for anything less 
than this. The three galaxies per region A and B is a conservative choice and is the 
least arbitrary limit that can be applied. Despite this choice, the signal is examined 
as a function of galaxy number in Supplementary Section 1. In sum: increasing 
the minimum number of galaxies in each region decreases the sample size and 
increases the significance of the measured signal.The final fiducial sample consists 
of 17,181 filaments, and 213,625 galaxies in total. 

A statistical significance to each measured value of Az,, must be obtained 
to ensure that the measured signal is not simply the result of stochasticity. In 
other words, we wish to ask the question “What is the chance that the observed 
distribution or Az,, could be reproduced given a random distribution of redshifts 
(in a given filament)?’ To answer this question, the following test is performed. 
The redshifts of all galaxies in a given filament are randomly shuffled, keeping 
the position fixed. This is performed 10,000 times per filament. In the poorest 


filaments with only 6 galaxies (say), there are less than 10,000 unique permutations. 


However, performing this test 10,000 times, although redundant, will not affect the 
determination of the significance. For each of these trials, the procedure described 
above is carried out: namely two regions are defined, their mean redshifts are 
measured and Az‘@?“°™ is computed. The observed Az, can then be compared 
with the distribution of random Azzn4™ on a filament-by-filament basis. The 
chance that a measured Az, is consistent with random can then be quantified by 
where (how many standard deviations from the mean) the measured Az,, is in the 
distribution of Az@ndo™, 
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